Kinetics of macroion coagulation induced by multivalent counterions 
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Due to the strong correlations between multivalent counterions condensed on a macroion, the net 
macroion charge changes sign at some critical counterion concentration. This effect is known as the 
charge inversion. Near this critical concentration the macroion net charge is small. Therefore, short 
range attractive forces between macroions dominate Coulomb repulsion and lead to their coagulation. 
The kinetics of macroion coagulation in this range of counterion concentrations is studied. We 
calculate the Coulomb barrier between two approaching like charged macroions at a given counterion 
concentration. Two different macroion shapes (spherical and rod-like) are considered. A new "self- 
regulated" regime of coagulation is found. As the size of aggregates increases, their charge and 
Coulomb barrier also grow and diminish the sticking probability of aggregates. This leads to a slow, 
logarithmic increase of the aggregate size with time. 

PACS numbers: 82.30.Nr, 87.15.Nn, 82.70.Dd, 87.14.Gg 



I. INTRODUCTION 

Water solutions of strongly charged particles 
(macroions) with multivalent counterions with large 
charge Z (Z-ions) are important in physics, chemistry, 
biology, chemical engineering and environmental science. 
Colloidal particles, charged membranes, double helix 
DNA, actin and other polyelectrolytes are examples 
of different macroions. Multivalent metallic ions, den- 
drimers, charged micelles, short DNA helices or other 
short PE can play the role of Z-ions. 

We concentrate here on strongly asymmetric solutions 
in which size and charge of macroions are much larger 
than those of Z-ions. As the simplest example, we con- 
sider macroions as negatively charged rigid spheres with 
charge —Q and radius R in solution with compact posi- 
tive Z-ions with the size a <^ R and charge 1 <^ Z <^ Q. 
This can be a solution, of positive latex particles with very 
short DNApjheliceaj'El or hematite particles with poly- 
acrylic acidcl. 

In such solutions, each sphere adsorbs many Z-ions. 
They strongly repel each other at the surface of the 
sphere and form a strongly correlated two-dimensional 
liquid reminding Wigner crystal. When a new Z-ion ap- 
proaches this liquid, it repels nearest Z-ions, creates a 
correlation hole or an oppositely charged irnage, which 
provides attraction of Z-ion to the surfaceo in addition 
to what mean field theories predict. Therefore, when con- 
centration of Z-ions, c, reaches some critical value co(s), 
(which depends on concentration of macroions particles 
s) the net charge of each macroion, Q*, (which includes 
all adsorbed Z-ions) flips its signQ. On Fig. 1 the "neu- 
trality" line co(s) is shown in plane (s, c) together with 
two signs of Q* which it separates. 

Correlations between Z-ions on the sphere surface also 
lead to another interesting effect. The correlation energy 
per Z-ion is lower at higher Z-ion concentration. This 
means when two spheres touch each other, at the place 
of contact where the Z-ion concentration doubles, cor- 




FIG. 1: Phase diagram of a solution of negative spheres 
with compact Z-ions in plane of their concentrations (s,c). 
The dotted line corresponds to the isoelectric composition 
c — sNi , where Ni is the number of Z-ions needed to neutral- 
ize one sphere. The dashed line corresponds to the concentra- 
tion of Z-ions, c = co(s), where net charge of a sphere with 
adsorbed Z-ions crosses zero. The two solid lines define the 
external boundary of domain where spheres coagulate. Plus 
and minus are the signs of the net charge of free spheres above 
and below the dashed curve. 



relation energy is gainedB This energy gain leads to an 
attraction between spheres which is stronger than van der 
Waals forces and makes possible their coagulation. Coag- 
ulation is a key part of many industrial processes such as 
paper production, extraction of minerals, proteins and 
other macroions from solutions, or treatment of waste 
waters. On the other hand, there are many other cases 
where coagulation should be avoided. Delivery of short 
modified DNA molecules DNA adsorbed on the surff 
of positively charged latex particles is a good exampL 

Because colloidal solutions are stable due to the 
Coulomb repulsion between particles, coagulation is usu- 
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ally achieved by a large concentration of monovalent salt, 
which screens out these charges. When the concentra- 
tion of salt grows beyond the coagulation threshold noth- 
ing new happens and macroions stay coagulated. If, in- 
stead of monovalent salt, we deal with Z-ions (add Z:l 
salt) then charge inversion changes this situation. Co- 
agulation happens only in the range of concentrations, 
which is close to the neutrality line (shown by gray in 
Fig. 1), so that the Coulomb repulsion is so weak that 
the free energy is gained when spheres touch each other. 
Thus when, say, the concentration of Z-ions, c, grows at 
fixed s, coagulation happens at some critical concentra- 
tion, Cc, while at a larger concentration, Cd, aggregates of 
spheres dissolve, because spheres acquire large positive 
net charge. In our previous paperd, we derived expres- 
sions for Cc(s) and Cd{s). These curves are plotted on 
Fig. |l|. Remarkably, phase diagrams of this type were 
discovered experimentally more than half a century ago 
for complexes of strongly oppositely charged proteinsQ. 
However, in Ref. ^ only the charge of large aggregates of 
macroions (coacervate droplets) was measured and the 
net charge of a single macroion was not discussed. No 
explanation was proposed for charge inversion of aggre- 
gates. In Ref. H we suggested an equilibrium theory of 
the phase diagram of Fig. |l| and showed that aggregates 
of spheres and isolated spheres change sign of their net 
charge at the same line co(s). 

In this paper we go beyond equilibrium statistical 
physics and study kinetics for different domains of the 
phase diagram. There are two main slow processes in the 
problem: charge inversion and coagulation. The first one 
is slow in the upper left area of the phase diagram when 
free macroions are strongly overcharged. In this case, the 
growth rate of the macroion positive charge is limited by 
a large repulsive Coulomb barrier for new Z-ions. Acti- 
vation above this barrier is necessary for a Z-ion to come 
close enough to macroion in order to feel attraction to its 
image in the strongly correlated liquid (SCL) of already 
adsorbed Z-ions. 

In the coagulation domain surrounding the neutral- 
ity line charges of macroions are relatively small so that 
charge inversion is a fast process. On the other hand, co- 
agulation of macroions can be very slow and take hours. 
Suppose we mix latex spheres with Z-ions in such con- 
centrations that the corresponding point is in the gray 
area of the phase diagram and watch how the mass of 
aggregates grows as a function of time. Although the net 
charge Q* in the gray area is smaller than the bare charge 
Q, its absolute value can be much larger than Z and the 
Coulomb barrier of repulsion of two macroions can be 
much larger than the barrier for charge inversion, which 
is created by the Coulomb interaction of a Z-ion with 
the macroion net charge Q* . Therefore, in this paper, 
we assume Z-ions are always in equilibrium and study 
the kinetics of macroion coagulation only. 

We first concentrate on the short time kinetics in which 
doublets of macroions appear. We calculate the Coulomb 
barrier between two approaching macroions as a function 



of the distance from the neutrality line. We measure this 
distance by the variable 
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where the critical concentration 



co(s) Co + sNi. 



(1) 



(2) 



Here Ni is the number of Z-ions needed to neutralize one 
sphere. The concentration cq is the concentration of free 
Z-ions which|-js in equilibrium with neutralized spheres. 
It is equal tcQ Cs exp(/i5Ci//cBT), where Cs is the three 
dimensional concentration of the surface layer of the SCL 
of Z-ions and ^scl is the chemical potential of Z-ions in 
SCL. It is negative and |/^scl| ^ ksT so that cq is very 
small. 

In Ref. H, the net charge of the sphere was shown to 
be proportional to U near the neutrality line: 



= CU 



(3) 



where C is the capacitance of the macroion coated by 
the "metallic film" of SCL of Z-ions. For a sphere C — 
eR, where e is dielectric constant of water. Because U 
depends on the correlation chemical potential ^scl of Z- 
ions in the SCL, we call it the correlation voltage. From 
Eq. (|l|), one sees that if c > co(s) spheres are overcharged 
(positive), while in the opposite case c > co(s) they are 
undercharged (negative). Detail derivation of Eq. (||) 
can be found in Sec. II of Ref. || (See Eq. (26) of this 
reference) . 

Our main result is that the height of the Coulomb bar- 
rier between two spheres is 



(4) 



where a is a numerical factor of the order of unity (at 
week screening, a — 0.3, at very strong screening, a = 1). 
Thus the rate v at which doublets appear has the form 



v{U) = iyoe^p{-aCU^/kBT) 



(5) 



Exponential factor of Eq. (^ is, of course, the probabil- 
ity of activation above the Coulomb barrier between two 
macroions. 

Eq. (^) suggests that the Coulomb barrier between 
aggregates increases as their size grows. This leads to the 
decreasing sticking probability between two aggregates as 
their size increases. This, in turn, leads to the slowing 
down of the growth of the aggregate size with time. We 
call this the "self-regulated" aggregation and show that in 
this case an aggregate size increases only as a logarithmic 
function of time 



7^(^) ~ 7^l in(</ti). 



(6) 



where the size and time constant TZi and ti will be given 
in Sec. III. To the best of our knowledge, this logarithmic 
behaviour is reported for the first time. 
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However, "self-regulated" aggregation cannot continue 
forever. When aggregates are so large that their time 
of activation above the Coulomb barrier is longer than 
the time for one macroion to desorb from an aggregate, 
one enters a Lifshitz-Slezov (LS) regime of coagulation 
where the aggregate size increases linearly with time with 
a very long time constant. In this regime, aggregates 
gains size by adsorbing free spheres which desorb from 
other aggregates. 

Screening can diminish the Coulomb barrier substan- 
tially. At strong screening, the LS regime may never 
be reached. Instead the height of Coulomb barrier satu- 
rates when the aggregate size reaches the screening length 
Vs. In this case, the "self-regulated" aggregation is fol- 
lowed by the reaction limited aggregation where Coulomb 
barrier is constant for any aggregates with size greater 
than Ts and the average aggregate size increases roughly 
quadratically in time. 

There have been various experiments studying the ki- 
netic of macroion coagulation induced by Z-ions, suck 
as latex particles complexed with short DNA segments ll 
or hematite particles complexed with polyacrylic acidE. 
In both experiments, the author observed an exponen- 
tial increase in the rate of coagulation when the neu- 
trality line is approached. At the same time, as the 
Z-ions concentration increases across this line, the elec- 
trophoretic mobility of the macroions changes sign sug- 
gesting a charge inversion effect. When Z-ions concentra- 
tion are much larger or smaller than the neutrality com- 
position, macroions are under- or overcharged and the 
Coulomb barrier exponentially diminishes their sticking 
probability. These results are in qualitative agreement 
with Eq. (||). However, instead of ln[z^(J7)] oc C/^, the 
authors of Ref. || observed a plateau in the aggregation 
rate around the neutrality point U = 0. This may be due 
to the limited time resolution of the experiment at the 
initial fast stage of coagulation. 

The paper is organized as following. In Sees. II and 
IV, we study the kinetic barrier for two approaching 
macroion with spherical shape and rod-like shape at a 
given Z-ion concentration. In section III, we discuss 
different stages of macroion coagulation, namely how 
the coagulation rate crossovers from a diffusion limited 
regime to the new "self-regulated" regime and finally to 
LS regime. In the conclusion, wc summarize our results. 



II. COULOMB POTENTIAL BARRIER 
BETWEEN TWO APPROACHING SPHERES 

Let us start by calculating the Coulomb barrier be- 
tween two approaching spheres. Because the charge of 
each sphere is not fixed but self-adjusts (by releasing or 
absorbing Z-ions) according to their positions, one has 
to calculate self-consistently the Coulomb repulsion be- 
tween spheres and their charges at a given separation. 
The capacitor charging picture of Eq. offers a very 
convenient way of doing this. Indeed, the voltage U de- 



pends only on the concentration of bulk Z-ions and the 
surface charge density of the spheres (through the chemi- 
cal potential hscl) and, therefore, is constant for a given 
c and sphere surface charge density a. Because the SCL 
of Z-ions on the macroion surface behaves as a metal, 
one can view the system of spheres and their aggregates 
as a system of conductors under a constant charging po- 
tential U . Thus, one can calculate the net charge of not 
only spheres but also their aggregate of any size by us- 
ing the appropriate capacitance Caggrogatc instead of C. 
In the same way, the kinetic Coulomb potential barrier 
between any two aggregates can also be calculated. 

Because, Z-ions are much more mobile than spheres, 
the process of charging up the spheres is much faster than 
coagulation and at any instance during the coagulation 
process, the Z-ions distribution is in equilibrium. With 
this assumption, one can write the Coulomb interaction 
energy between two spheres when they are at distance r 
from each other as: 

where Cii(r), C22(?') and Ci2(r) are, respectively, the 
self capacitances of the spheres 1 and 2 and their mutual 
capacitance. The capacitance of an isolated sphere of is 
eR. 

When dealing with a system under constant charging 
potential, the free energy F{r) of the system must include 
the work of the source (in our case, the population of free 
Z-ions) to maintain this potential 

2 

F{r)^V{r)-Y,C,,{r)U^^-V(r). (8) 

ij' = l 

It is the lowering in the free energy compared to the 
reference system of two completely neutralized spheres 
with the rest of Z-ions free. 

The capacitance Cn, C22 fl-nd C12 has been calculated 
for two spherical conductorsEl: 

00 

Cii(r) = C22(r) = £i?sinh/3^sinh"^[(2n- l)/3], (9) 

n=l 
00 

Ci2(r) = -£i?sinh/3^sinh"\2n/3), (10) 

n=l 

where /3 > satisfies cosh (3 — zl^R- 

Substituting Eqs. (^) and ( p^ into Eq. (^, one can 
easily find the cost in the free energy —V{r) in moving 
two spheres from infinity to the distance r. The result is 
plotted by the solid line in Fig. ||. For comparison, the 
potential barrier for the case when the spheres keep their 
charge q* — C{oo)U — eRU fixed when approaching each 
others is also plotted. The maximum of the potential 
barrier is located very close to the distance r — 2R where 
two spheres touch each other and is equal to 

V^acc ^ O.SeRU^ (11) 
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FIG. 2: The Coulomb potential barrier between two spheres 
(the solid line). For comparison, the Coulomb potential bar- 
rier for two spheres with fixed charge is plotted (the dashed 
line). Note that, although not plotted, at distances very close 
to r /2R = 1, the interaction energy drops below zero due to 
the short range correlation attraction. 



As one can see, at the maximum of the potential bar- 
rier, the self adjustment of the sphere charge reduces this 
barrier height by about 40% compared to the barrier of 
Q.heRU^ one would get if the sphere charges were kept 
fixed upon approaching. 

Due to this Coulomb barrier, the expression for the 
rate at which spheres come to each other to form doublets 
contains an exponential factor related to the probability 
of activation above this barrier as shown in Eq. (^). 

In the above calculation, it is assumed the Coulomb 
potential is unscreened. In reality, there is always a finite 
concentration of monovalent salt in water solution, which 
leads to the screening of Coulomb interaction at distance 
larger than the Debye-Hiickel screening length . If Vg 
is larger than the averaged distance between Z-ions at 
the surface of an aggregate, the correlation between Z- 
ions remains unscreened. Therefore, the concentration 
Co, and correspondingly, the charging potential U remain 
constant. Thus, in this regime, screening influences the 
process of coagulation only through the change in the 
capacitance of each sphere. 

When the radius of the aggregate is smaller than the 
screening radius, R < r^, the Coulomb interaction is not 
screened when the aggregates touch each other, therefore 
the above result remains valid. In the opposite case, when 
R > rg, the capacitance of each isolated aggregate is that 
of a plane capacitor with thickness and area 47ri?^, 
Ci.2(oo) = sR'^/ts. Using this planar capacitor approx- 
imation, it is not difficult to estimate the change in the 
total capacitance of the system when the two aggregates 
touch each other. Indeed, simple geometric calculation 
gives the reduction in the area of the planar capacitor 
when the two spherical aggregates touch each other is 
iirRrs. Thus the change in the capacitance of the sys- 



tem is AC = Cii -I- 2Ci2 + C22 - Ci(cx)) - C2(cx)) = eR. 
Thus, in this case the maximum of the potential barrier 
is 



Vrr, 



eRU^ 



{R»rs 



(12) 



Comparing this results with Eq. (11), one sees that, 
beside a numerical factor, screening does not affect the 
height of the potential barrier between two approaching 
spheres as long as Ts is larger than the average distance 
between Z-ions on the spheres. 



III. LATER STAGES OF COAGULATION 

In previous sections, we discussed the kinetic barrier 
when two macroions approach each other to form a dou- 
blet. This fast process happens at the initial stage of 
coagulation where the concentration of free macroions is 
large. Doublets coagulate with free macroions or other 
doublets to form bigger macroions aggregates. This sec- 
tion deals with later stages of coagulation. 

The problem of coagulation is generally complicated 
and can be solved exactly only in special cases. However, 
the main physics can be captured if one works with the 
dominant aggregate size Tl{t) (with concentration N{t)) 
at a given time assuming this typical aggregates carry all 
the mass of the macroions. 

Let us start from the initial stage of coagulation when 
the aggregate size are small, so that at small enough 
U (close to the neutrality line) the Coulomb barrier be- 
tween aggregates is smaller than ksT. In this case, the 
Coulomb barrier of Eq. ( pdj ) has little effect and the 
probability of sticking for two aggregates is of the order 
of unity. This is the regime of the well-known diffusion 
limited aggregation. The concentrationxtf the typical size 
aggregates decreases as inverse of timdlil: 



Nit) 



——^ ~ for t > Tdiff , (13) 

1 + Tdiff t 



where Tdifi rj/ksTs, s is the initial concentration of free 
macroions and rj is the viscosity of water. Thus, the rate 
of coagulation in this early stage is relatively constant 
independent of U and equal to the rate of coagulation 
at the neutrality line where the charging potential U is 
exactly zero and there is no Coulomb barrier. 

The regime of diffusion limited aggregation stops when 
the typical size of the macroion aggregates reaches such a 
size that the Coulomb barrier is larger than fc^T. Using 
Eq. (pd]), one sees that this regime is reached when the 
typical size equal 



(14) 



Assuming an aggregate has a fractal dimension df (in 
practical situation, ~ 2) so that the number of spheres 
in it is {TZ/R)'^f , the time at which TZi is reached is 



ti = Tdiff (fcsT/eC/2i?)" 



(15) 
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When t > ti, the Coulomb barrier Vmax — sTZU'^ has 
to be taken into account. In this regime, the dominant 
contribution to the rate of coagulation comes from the 
tunneling of aggregates through the Coulomb barrier. 
The change in the_concentration of typical aggregates 
obeys the equatioiO: 



dN{t) 
dt 



^ -KDTZ^N^exp{~V,nax/kBT) (16) 



where the numerical factor has been dropped on the 
right hand side and the factor k satisfies = 
~{d'^V{r)/dr'^)max/'2kBT. Since the typical width of the 
Coulomb barrier is TZ, k ~ y/Vmax/Tl^kBT. Using the 
Stoke formula, D ~ kBT/rflZ, one can solve Eq. ( [iq ) as- 
suming that the typical size aggregates consume all the 
mass of the macroions: 



7^ = 7^l(A^l/iv)l/''^ 



(17) 



where iVi = N{ti). The concentration of typical size 
aggregates in this regime decreases logarithmically with 
time: 



N{t) ~ A^i 



In 



t - ti 
tidf 



(18) 



where Ni = N{ti). Eq. (|) for the size of the typical 
aggregates can be obtained easily using the relationship 

As of our knowledge, this slow logarithmic kinetics was 
never reported in literature. It is the result of the increase 
in Coulomb barrier when aggregate size grows. This di- 
minishes their sticking probability, which in turn slows 
down the kinetics from the linear size increase to a log- 
arithmic one. We thus call this regime "self-regulated" 
aggregation. 

In the regimes we considered so far, the increase in the 
aggregate mass is caused by a collision between two ag- 
gregates of size TZ{t). When the typical aggregate size 
becomes very large, the Coulomb barrier between two 
approaching aggregates becomes so high that the cor- 
responding coagulation rate (which is proportional to 
e:iqy{—Vmax/kBT)) becomes very small, a different and 
faster aggregate growth mechanism, namely the Lifshitz- 
Slezov (LS) one, comes into play. In this mechanism, 
the large aggregates do not collide into each other and 
the primary mechanism for aggregate growth is no longer 
due to real space diffusion. Instead, the change in their 
size is provided by the releasing and adsorbing of spheres 
which have much smaller kinetic barriers. This mecha- 
nism provides a faster coagulation process because the 
binding energy of a sphere in a very large aggregate and 
the Coulomb barrier between a very large aggregate and 
a sphere are finite and do not depend on the aggregate 
size. One|-ean think of LS mechanism as a diffusion in 
size spacalJ. Interestingly, the concentration of typical 
size aggregates in this mechanism also decreases as in- 
verse of time 



n 




' R 




IJ 



FIG. 3: The typical aggregate sizes TZx and 7?,2 at which ag- 
gregate growth mechanism crosses over from diffusion limited 
regime to "self-regulated" regime and then to Lifshitz-Slezov 
regime. 



where the time constant tls is much longer than Tdig 
because it is proportional to the concentration of free 
macroions in solution 



(20) 



Here E is the binding energy of a sphere in an aggregate 
which does not depend on the aggregate size. 

One can easily find the size of the typical aggregate 
at which LS mechanism is important by comparing the 
detaching time t^s of a small macroion from the ag- 
gregates with the activation time for two aggregates to 
go through the Coulomb barrier which is of the order 
exp{-enU^/kBT). This gives 



eUU^ 



E. 



(21) 



N{t) cx TLsA, 



(19) 



Thus the typical size at which LS regime starts is: 

7^2 = E/eU^. (22) 

A phase diagram of different regimes of coagulation is 
shown in Fig. |[ With increasing time, TZ grows along 
a vertical line starting from TZ — R. Below the dotted 
line which represents TZi{U), the coagulation is diffusion 
limited and the aggregate mass increases with time lin- 
early as i/rdiff. Above the solid line which represents 7^2, 
the coagulation process is of LS nature and the aggregate 
mass also increases linearly with time ^/tls- In between 
these two lines is regime of "self-regulated" aggregation. 
In this regime, the aggregate size increases logarithmi- 
cally. 

Above, we dealt with the unscreened Coulomb poten- 
tial between two aggregates. This is valid if the screening 
length of the solution is very large (r^ > 7^2) ■ For smaller 
screening length, when the typical aggregate size reaches 
Ts, due to the fractal Ujature of aggregates, the Coulomb 
barrier stops increasingpil. In this case, LS regime cannot 
be reached. One goes from the "self regulated" regime to 
the regime of reaction limited aggregation with constant 
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FIG. 4: Sketch of the dependence of the size TZ of typi- 
cal aggregates as function of time t. The sohd Une corre- 
sponds to the case of very strong screening, Ts < TZi, where 
Coulomb barrier is smaller than ksT and aggregation is dif- 
fusion limited (TZ oc yJTprdm)- The dashed line corresponds 
to intermediate screening, TZ\ < Ts < TZ2, where one goes 
from diffusion limited aggregation to self-regulated aggrega- 
tion {TZ oc \n.{t/ti)) and to irreversible reaction limited aggre- 
gation {TZ oc i/rrcact)- 



Coulomb barrier and almost zero probability of macroion 
detaching from the aggregates. One can calculate the 
change in the concentration N in this regime using the 
Eq. (0). 

Denote by ts the time at which the typical aggregate 
size reaches Vg (from Eq. (^), ts ~ tiex.p{srsU'^ /kBT)). 
At t > ts, the exponential factor in Eq. ( p^ becomes 
constant, exp{—ersU^/kBT), and so does the factor k. 
The solution of this equation shows a fast decrease in the 
concentration of typical aggregates with time: 



Nit) = N{ts 



t — df — 1 



'^rcact 



df/a-df) 



(23) 



where Tpoact is an exponentially long time constant: 



eMersUykBT) 



kBTN{ts) ^ersUykBT 

In a typical situation, 2± 2, the number of spheres in 
the typical aggregates (which is inversely proportional to 
N) increases as t^. Although this is a very fast kinetic, it 
is still slower than an j expo nential growth of the cluster 
size suggested recentl5it2lEJ. We currently do not have a 
clear understanding of the origin of this theoretical result. 

At even smaller screening length, < TZi, the 
Coulomb barrier between two aggregates never becomes 
larger than fc^T and one always stays in the regime of 
diffusion limited aggregation. 

The evolution of the typical aggregate size as a function 
of time is plotted in Fig. |4|for Vs < Tl2- 



a) 




b) 



FIG. 5: Different paths of coagulation for two rigid rods, (a) 
The path studied in Ref. ^ (b) The path which according to 
our paper has lower kinetic barrier. 



IV. AGGREGATIONS OF ROD-LIKE 
POLYMERS 



In this section, we would like to discuss the kinetic bar- 
rier for the coagulation process when the role of macroion 
is played by a rigid polyelectrolytc (PE), such as DNA. 
We assume the length of PE molecules is smaller than 
its persistence length so that each can be considered as 
rigid rod. Due to the anisotropy of this problem, there is 
a number of different paths of aggregation of rods. 

The authors of Ref. || studied the kinetic barrier be- 
tween two approaching rods as a function of their orien- 
tation when the rods approach each other in the direction 
perpendicular to their bodies (Fig. ||a). We believe the 
kinetic barrier associated with this way of approaching 
is too large because of their large electrostatic repulsion 
when placed side by side, especially when the screening 
of solution is week. In this section, we propose another 
path of approaching of two rigid rods which has a lower 
kinetic barrier. Namely, the two rods approach in such 
a way that their centers of mass lie on a line parallel to 
their longer axes (See Fig. ^3). One can easily see that, 
in this way, charges of the rods are kept further from each 
other than in Fig. ^ja. Thus for the paths of Fig. the 
kinetic barrier between the rods is lower. 

We can calculate the Coulomb barrier between two ap- 
proaching rods as in Sec. II, starting from the idea that 
both rods are under constant charging voltage U. The 
capacitance of a metallic rods of length L and radius 
a is C = eL/2ln{L/a) for a weakly screening solution 
{rs >> L) and is C = eL/2ln{rs/a) for stronger screen- 
ing [L » Ts » a). Thus the change in the capacitance 
of the system when the two rods ends start touching other 
is AC ~ eL/\n^{L/a) for the weak screening case, and 
is AC ~ erg/ In^irs/a) for stronger screening. 

As the rods start to overlap, the short range attraction 
between the rods starts to come into play, which partially 
compensate the increases in Coulomb repulsion and re- 
duced the kinetic barrier. If the length of the overlap 
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segment is x and if the attraction energy is — e per unit 
length, the free energy of the system can be written as 
-AC(a;)C/2 - xe, which is x^eLU"^ / \n{L / a) - e) for the 
weak screening case and is x(eXf/^/ ln(rs/a) — e) for the 
strong screening case. Because, in the final stage when 
the two rods lie parallel to each other, the attraction 
energy is assumed to win over their Coulomb repulsion, 
e > eLU"^ / \n{rs/a), one easily sees that when the rods 
start to overlap, total energy start to decrease. Thus the 
maximum of the potential barrier is at about the distance 
at which the rods start to overlap. 

„ _ r eLUyin^{L/a) for r, > L , . 

Vraax - <y {/2/ ^^^2 ^ ^ {^V 

Obviously, this potential barrier is much smaller than 
that in the case the rods approach each other in the di- 
rection perpendicular to their length. In the latter case, 
the maximum of the potential barrier is at the distant 
where they touch each other side by side and is equal 
to eLU'^/\n{L/a) for > L, and eLU^ / \n{rs/a) for 
rs < L. 

At later stages of aggregation, coUinear approaches 
should dominate as well. This can explain why in many 
real and numerical experiments, very elongated struc- 
tures are seen. 



V. CONCLUSION 

In this paper, we provide a general framework for cal- 
culation of the Coulomb potential between two macroions 



in solution of multivalent counterions and study various 
stages in the coagulation of macroions. The capacitance 
interpretation of Eq. (||) proves to be very useful in calcu- 
lating the net charge of any macroions and aggregates of 
any shape. It also helps to easily calculate the Coulomb 
barrier between two approaching macroions. Using this 
equation, we are able to calculate and compare different 
paths of coagulation for rod-like macroions. 

We discuss several stages of coagulation in time. In 
low salt, the coagulation process goes from a diffusion 
limited regime to a "self-regulated" regime and finally to 
the regime Lifshitz-Slezov kinetic. In the "self-regulated" 
regime, as the aggregate size increases, their Coulomb 
barrier increases diminishing their sticking probability 
and slowing down the kinetic. As a result, the aggre- 
gate size increases as a slow logarithmic function of time 
instead of the standard linear relationship. At higher salt 
concentration. Coulomb barrier is screened and stop in- 
creasing after the aggregates reach a certain size (of the 
order of the screening length rs). In this case, Lifshitz- 
Slezov regime cannot be reached. Instead, one reaches 
a reaction limited regime where the Coulomb barrier is 
constant and the aggregate size increases quadratically in 
time. At very high salt concentration, one even cannot 
reach the reaction limited aggregation regime and always 
stays in the regime of diffusion limited aggregations. 
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